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Abstract 

Thanks to averaging processes and F-convergence techniques, we 
are able to link a microscopic description of ferromagnetic materials 
based on spin lattices and their mesoscopic description in the static 
framework for the three fundamental contributions: exchange, magne- 
tostatic and external field. The results are in accordance with the clas- 
sical continuous description of ferromagnetic phenomena and justifies 
it. This work is a seed towards a dynamic description of ferromagnetic 
materials. 

1 Introduction 

The continuous description of ferromagnetic materials has been introduced 
since the 60's via the micromagnetism model developed by W.-F. Brown 
This model, based on a thermodynamical description of ferromagnetic phe- 
nomena, has proved its efficiency in numerous works via relevant simulations 
( [21 El El H] ) • Nevertheless, several problems persist in the description, from 
the thermic effects to the magnetostrictive behaviors. These problems are 
very sharp and, in order to understand their modeling, one need to under- 
stand the microscopic behavior of magnetization (atomic scale) and the link 
between this scale and the mesoscopic scale (continuous magnetic matter 
scale of ferromagnetic effects). In this paper, we focus on the beginning of 
this program: the link between a microscopic description of ferromagnetic 
materials and their mesoscopic description in the static framework for the 
three fundamental contributions: exchange, magnetostatic and external field 
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[HITllH]. Here, the microscopic scale designates the atomic scale where atom 
nuclei are assumed to be point electric charges bearing one magnetic moment 
induced by the atom's electronic cloud; the mesoscopic scale designates the 
continuous description of matter for which the ferromagnetic effects are sig- 
nificative. 

At the microscopic scale, we describe the material gular periodic 

spin lattice intersected with the magnetic domain (a bounded open set of M^). 
Section [2] is dedicated to the mathematical description of the microscopic 
model (the spin lattice) and to an averaging process towards a mesoscopic 
models which leads to a constant norm magnetic field, in accordance with 
usual models of micromagnetism. The microscopic energies are introduced 
and several modeling hypotheses are set. The main hypothesis is induced by 
the adiabatic behavior of the Heisenberg energy compared to the global ferro- 
magnetic energy. This hypothesis gives a constraint on neighboring magnetic 
moments. In fact, this constraint is verified by a set of minimizers for a given 
lattice. 

The energy induced by the Heisenberg interaction is more difficult to 
treat. Section [3] addresses the study of this contribution and, in particu- 
lar, its asymptotic behavior for sequences of lattices verifying the modeling 
hypothesis. This hypothesis ensures compactness which allows to use F- 
convergence tools in H^. The limiting energy constructed from the discrete 
magnetization is the exchange energy (Theorem [1]). 

In Section IU we introduce the demagnetization energy both for discrete 
lattices and the continuous model and finally obtain a convergence result for 
the sum of the Heisenberg and demagnetization contributions (Theorem [2]). 

2 Mathematical descriptions of a spin lattice 
2.1 Atomic lattice description 

We consider a collection of spins which are located on the nodes of a periodic 
lattice C in the M.'^ space {d = 1,2,3) with mesh size a > 0. In the scope 
of this paper we will restrict to the case of ID, square or cubic lattices, £ is 
simply aZ*^, but we can think of more complex lattices. Here all the nodes 
play the same role to ensure a unique definition of neighbors. In the opposite 
case a multi-species model should be used. 
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The nodes are indexed by z e N and we denote by Xi the zth spin location 
and the corresponding spin value (magnetic moment). The norm of these 
magnetic moments are scaled to the unit value and therefore for all i e N, 
/Xxi e S^, where is the unit sphere of M^. 

Instead of describing a collection of magnetic moments, we can gather all 
the values in one single vector field /i defined by 

where is the Dirac delta function centered at Xi. 
2.2 Scaling 

We want to obtain an homogenized model of the spin lattice, i.e. give a 
description when this lattice is seen from far. Instead of really doing this, 
we will perform some dual transformation, i.e. consider only nodes that are 
included in some fixed bounded domain fi, and shrink the lattice (as shown 
on Figured] for d = 2). More precisely, we suppose that e and for all 
n e N*, using the homothety hn{x) := x/n, Vx e M'^, we define 

• Cn = hn{C), the shrunk lattice; 

• j(^n,Q = Cn <~^^, the nodes of the shrunk lattice that belong to ^2; 

• fin ^ (S^)^"'", the shrunk vector field. 
We notice that for all y eW'-, 

where fin,x = l^h-^{x)- assume that f2 has a sufficiently regular boundary 
in order that the number of nodes belonging to Cn,n is 

i^Ln,n = Cn\l + 0{^)). 

where C is a constant which only depends on a and Vt (which are constants 
of our problem). 
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(a) (b) (c) 

Figure 1: Scaling of a 2D spin lattice: each sub-plot represents the domain 
and (a) the square lattice £; (b) the shrunk square lattice £2; (c) the 
homogeneized lattice in O as n ^ co. 



2.3 Regularity assumptions 

In order to pass to the limit as n ^ 00, we have to assume that the magnetic 
moments are locally almost aligned. The definition of locality is given by an 
integer multiple A; e N* of the shrunk mesh size ajn. 

We define a first regularity assumption that only depends on the distance. 
For all X G R'^ and r > 0, we denote by B{x, r) the ball of center x and radius 
r in R*^. 

Hypothesis 1. For all n eW , there exists Cn > such that 

( ka\ 

where Cn = 0{l/n^). 

We are indeed interested in the limit when we have a more and more dense 
lattice of more and more aligned spins. We therefore perform a diagonal 
process and correlate n and the spin alignment. 

To define the averaging process we will also need to assume the Q has the 
uniform cone property. 

Hypothesis 2. There exists an angle a and a radius r, such that for ally e fl 
there exists a cone Cy of angle a with vertex at y such that B{y,r) nCy a fl. 
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2.4 Partitions adapted to the lattices 



Let us first work on the initial lattice. To this aim, we define a partition of 
unity {px)xekc adapted to the dilated lattice kC. Since our lattice is uniform 
and all the nodes are equivalent, all the px are equal up to a translation (see 
Figure [2]), i.e. there exists p* e C^f (M'^;M) such that 

VxG A;/:, VyeM^ px{y) = p\y - x). 




Figure 2: Partition {px)xekc in dimension 1. 

Without lack of generality, we assume that p* > and supp p* g B{ka). 
By definition of a partition 

vi/gm^ Y,p-^y)= 2]p*(i/-^) = i- 

xekC xekC 

Let rik be the number of nodes in B{0, ka), which is e.g. rik = {2k + 1)'^ for a 
cubic lattice. If we now sum over all the x e C,we have the same translation 
property and 

Vz/gM^ Y,P*{y-x)=nk. 

xeC 

Defining p = W^^p*, we have 



Vt/GM^ }]p(y-x) = l, 



xeC 



and therefore a continuous partition of unity on 

Now, we adapt this partition to the shrunk lattice Defining 

Vn G N*, pn{x) = piK^x), 
we have a continuous partition adapted to the shrunk lattice (see Figure [3]) 

yy G R'^, 2j Pn{y -x) = 1 and supp p„ cz B{0, — ). 

xeCn 



ajn 



ka/n 



Figure 3: Partition {px)xec„ in dimension d = 1. 
Since Vp* is uniformly bounded (i.e. 0(1)), then Vp„ is uniformly 0{n). 

To define an averaging process we will have to use a truncated partition 
of unity, namely 

^n{y) = 2 pn{y-x). 

We clearly have ^ $„(y) ^ 1 and as a finite sum of functions, e 
C"^(M'^;M). We also have a stronger result, namely $«(?/) is bounded from 
below uniformly in n and y e Q: there exists 6 > and uq e N* such that 

Vyefi, \fn^no, 6 ^ $„(?/) ^ 1. 

This result stems from the "cone property" (Hypothesis [2]). Besides is 
uniformly 0(n). 

2.5 Averaging process 

From the sequence {fj.n)nem of vector fields, we now define a new sequence 
{'nin)n^no '^^ regular vector fields on VL: for all n ^ rio, we define 

1 

Remark 1. Since pn and e C'^(M'^;M), and $„ is hounded from below by 
b > 0, it immediately follows that rUn e C°°(f2; M^). In particular, since Q is 
bounded, rUn e H^{Vl;B?). 

Lemma 1. The sequence {mn)n-itno is bounded in L°°(f2;]R^) and L^(f2;]R^). 
Proof. Recall 
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hence, for all y e f2, 



We clearly have 



< 



^ S K.|Pn(y-x) = |^ = l, 



||"^n||L«(Q;R3) ^ 1, 



□ 



Proposition 1. Under Hypothesis [21 {mn)n^no is a bounded sequence in 



Proof. For all y e we write 



V$.(|/) 



50 m, 



Let e n supp p„(?/ — •), we can write 



Hence 



X ^ Pn{V - X) ® Pn,xy 
+ 2 - a;) ® {pn,x - Pn,xy) 

= V^niy)®Pn,Xy 

V$„(?/) , , V$„(|/) 
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x) 



J 



+ 



+ 




■y 



In both sums there is only an 0(1) number of terms for which p„ or Vp„ 
is non zero. In the first sum, by Hypothesis [H \iJLn,x — Pn,xy\ = 0{(n'^) = 
0{l/n). Besides |V$„(?/)| = 0{n) and is bounded from below. The 
first term is therefore 0(1). In the second sum \fin,x — fJ'n,xy\ = 0{l/n) and 
^ Pn{y — x) = 0{n). Hence there exists O' > such that, for dA\ y e fl and 
n ^ no, |Vm„(?/)| ^ C 

Finally, Q being bounded, hence {\\'Vmn\\L2(n-K3d-^)n^no is a bounded se- 
quence. □ 

2.6 Asymptotics towards a mesoscopic model 

Since H^(Q;M.^) is weakly compact. Proposition [T] implies that there exists 
m G ifl(^];M3) such that 



From now on, we have to assume that Q is compact (closed) and has a 
piecewise boundary, to ensure that this implies that this convergence is 
strong in L^(f2; M^). 

Proposition 2. Under Hypothesis [21 m has a constant norm equal to 1 
almost everywhere on Vt. 

Remark 2. We recover here a property of the magnetization field in Brown's 
model of micromagnetism J4^, where the constant norm is assumed. 

Proof. We have seen that for all n ^ and y e Q, |m„(?/)| ^ 1. We want to 
show that it is also bounded from below and pass to the limit. For all y e Q, 



rrin m weakly in H^(Q] M^). 
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since \Hn,x\'^ = 1, 

\^n{y)\'^ = Yj I^n,xpn{y-X) 

= -^rj-yl Yj tJ-n,x ■ IJ^n^x' Pn{y - X)pn{y - x'). 

Since supp p„ cz B{0, ka/n), the sum runs indeed on Cn,n n B{y, ka/n): 

\mn{y)\^ = ^] \2 I] Pn,x- Pn,x' Pn{y - X)pn{y - X'). 

Hypothesis [T] imphes that 

IxiCL 

n 

and we also have that 

VyeJ^, J] pn{y - x)pn{y - x') = ^n{yf, 

x,x'eC„n 



therefore 



hn(?/)|' ^1- 



In the L (fi; R ) hmit, we therefore have \m\ = xn almost everywhere. □ 



3 Towards the exchange energy 
3.1 Heisenberg interaction 

The interaction of the spins is described by the Heisenberg energy, which can 
be written as follows: 

^x,y '^-^x,y{Px ' Py l); 

where A^^y ^ only depends on the distance between x and y. We make the 
assumption that each node x only interacts with its neighbors A^^;. 

Let Vo be the set of neighboring nodes in interaction with the node 
(0, 0, 0) via the Heisenberg energy. In the shrunk lattice, we will restrict 
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the computation of this energy to elements x and y m. e Cn^n- Since the 
lattice is homogeneous, the neighbors of any given node x e Cn,n can be 
deduced from the definition of A^o^ Nn,n,x = {{x} + /i^^(A^o)) ^- We also 
assume that (a) Nq (and hence Nn,n,x) is a finite set; (b) A^o -6(0, ka) (i.e. 
Nn,n,x <= B{x, ka/n)). 

We therefore define the node energy by 

&x = Gx,y ~ ~2 S ^x,y{fJ'X ' /^y ~ !)• 

In what follows we will restrict the study to 3D cubic lattices. Dimensions 
1 and 2 are also possible to treat in the same way. The fact that the lattice 
is cubic is used in the explicit computations, but our proof may be extended 
to other classes of regular lattices. We also consider as neighbors only the 6 
closest ones on the cubic lattice, which at at the same distance. Since A^^y 
only depends on the distance, which is now equal for all the neighbors, we 
can set 

'A>0, ifyeN,; 
0, else. 



Ax. 



y 



For all n e N*, we define the exchange energy of the domain Q associated 
to the spin distribution /i„ summing up the node energies scales to 
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3.2 Spaces and convergence 

We define the space sequence {Wn)nEm by Wn = (M^)^"'". In the sequel, we 
set W = n„eN* 

Hypothesis 3. Let fi = (/i„)„eN* ^ IlneN* There exists c > such that 



c 



Vn G N*, \/x e Cn,n, Vy G Nn,n,x, W^x - /^n,y| ^ 
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except possibly for a subset ln,n cz Cn,n such that i^{ln,n) = 0{n), and there 
exists a sequence {cn)n^i e such that Cn > and 

Vn G N*, Va; G ln,n, G Nn,n,x, \^i■n,x - fJ'n,y\'^ ^ Cn. 

More precisely, given the constant c, x e ln,n if there exists y e Nnn,x such 
that \^n,x — fJ^n,y\'^ > cjr? . Hypothesis [3] is a little weaker than Hypothesis [H 
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We denote by cz W the set of /x = (/i„)„gN* for which Hypothesis [3] 
holds and cz the set of ^t e for which \fin,x\ = !> Vn e N*, Vx e 

We want to define the convergence of elements of W towards elements of 
the limit space H^(Q;M.^). To this aim we define a partition of the domain 
c in tetrahedra, by groups of 5 tetrahedra (see Figure H] and Section 




Figure 4: Partition of a single 3D lattice cell. These five tetrahedra are ele- 
ments Tn, and T„ for the corner tetrahedra, and T* for the center tetrahedron 
(see Section [3.3.ip . 



Definition 1. We define the projection Pn '■ Wn — > H^(Q;M.^) such that 
PnifJ-n) is equal to on the lattice nodes and linear on each element of the 
partition. 

Definition 2. We say that fi e W converges to m e H^{'[l;'R^) if Pnil^n) 
converges weakly in H^(Q; R^) to m as n ^ cc. We denote this fin ^ ^■ 

Definition 3. We define the projection pn '■ C{Q; M^) — > Wn, such that for 
all m G C{VL\ R^), Pnijn) e Wn and for all x e Cn^n, ipn{i^))x = rn^x). 

3.3 Main result 

The exchange energy £„,ex is a functional defined on Wn- The main result of 
this paper is the following 

Theorem 1. Let {Hn)nem ^ . In the sense of the topology defined by 
Definition [H 

^n,ex ^ ^oo,exi 

n— >oo 

where 

m I — > 2A |Vm(x)|^dx. 
Jn 
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Moreover, this convergence is compatible with the unit norm constraint. 

Remark 3. The T-convergence result is two-fold fBjj: 

construction: for all m e H^{Q; M^), there exists i-i e W such that ^in ^ 
and limsup„^o3£„,cx(^n) ^ ^^oo,cx(m). 

lower semi-continuity: for all /x e W , such that \^n\ ^ B, fin ^ "^oo ^ 
H'^{Q;R^) an(i liminf„^oo^n,ex(^n) ^ ^oo,cx("^oo)- 

Proof. The proof of Theorem [1] sphts into many steps to which various lem- 
mas are devoted, the technical proofs of which are postponed. 

Lemma 2. Let (/i„)„gN* e , then there exists moo ^ H^(Q;M.^) such that 

n— >oo 

Now this function moo is only in i^f^(fi;R^), and not continuous, and 
we need to use pointwise values of this function. Therefore, for all e, we 
approximate moo by a C^{VL\ M?) which coincides with moo on a smaller domain 
fig. We denote by A the Lebesgue measure on R'^. 

Lemma 3. Let u e H^{VL\M.^), for all e: > 0, there exists an open set Uu,e 
and a function e C^(^7;R^) such that for all x e VL\uju,e-, Ue{x) = u{x), 

\{(jJu,e) < O'lT'd 

lim — < e. 

Corollary 1. Let u e H-^iQjM."^), for all e > 0, there exists an open set 
Qu,e ^ with piecewise boundary such that X{Q\Qu,e) ^ £ o^^^^ 

eCi(i7„,,;R3). 

Lemma 4. Let u e C^(n;R^), then {pn{u))nem ^ , Pn{u) > u and 

?i— >oo 

lim^^oo ^n,ex(Pn('ii)) = ^^oo,cx ("^i) ■ Morcovcr if u has a constant unit norm, 
then {pn{u))nem e W^. 

For all £ > 0, we can therefore associate to moo an open set fi^ on which 
it is and define the approximate energies: 

^n,eJMn) = - 2 2 ^l/^n.y " /^n.x forallnGN*, 

^^,cx("^oo) = 2A\ \Wm^{x)\^dx. 
Comparing S^exil^n) and £^^^{pn{moo j), we first show that 
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Lemma 5. For all e > 0, liminf„^oo ^^,ex(Mn) ^ ^^,^,ex("^oo)- 

Estimating the remainder of the sums and integrals on u^, we can then 
prove 

Lemma 6. hminf„^oo ^n,ex(^n) ^ ^^oo,ex(?^oo)- 



3.3.1 Proof of Lemma [2t Limit of a lattice of spins 

Lemma [2] is proved using an explicit computation of the projection P„(/i„). 
To this aim, we define, for all n eN*, 

• Tni the set of tetrahedra which form the partition of Cn.n', 

• Tn, the set of corner tetrahedra (4 for each mesh of the lattice, see 
Figure H]) ; 

• T*, the set of center tetrahedra (1 for each mesh of the lattice, see 
Figure H]) ; 

• En, the set of edges of the mesh of Cn,n (given by couples of the indices 
of the lattice nodes); 

• C„, the set of edges of elements of T*; 

• Sn, the set of outer surfaces of T„ (triplets {i,j,k) where e C„, 
{i,k) and {j,k) e En). 



Remark 4. Generically fi„ = yjreTn'^ £ ^- suppose that the geometry of 
Vt is such that X{VL\fLn) = 0{l/n), i.e. there are order tetrahedra covering 
the difference set. We set Pn{l^n) to zero on Q\Qn ■ 



□ 




'n 



Figure 5: Different sets of tetraedra and edges. 
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Step 1. Estimate of |P„(//„)(a;)pdx. 

Let X e Qn, then there exists t e Tn and x e r. Since Pn{iin) is hnear 
on T, then Pn{fin)ix) is a pondered mean of the ^n,x^ ^ i = 1, • • • ,4, where 
(xl ) e {J^n,ny are the vertices of r. Therefore |P„(/x„)(a:)| ^ 1. 

Hence 

\ \Pn{f^n){x)\''dx = f |P„(//„)(x)|Mx= 2 f l^nW(a;)|Mx 



2 A(r) = A(n„). 



Step 2. Explicit computation of |VP„(//n)(a^)pda;. 

We compute \W Pn{^n){x)\^dx exphcitly using the tetrahedron decom- 
position of the lattice. On each tetrahedron Pnil^'n) is hnear, and therefore 
its gradient is constant. 

If T G r„, we can construct an orthogonal system using the lattice nodes 
{xl,xl,xl,xl) e {Cn,ny, for example (fin^xj - ^n,a;[)i=2,...,4- Each direction 
yields one component of the gradient. Since the length of the edges are a/n, 
the component of the gradient are the {nn,xr — lJ'n,xi)n/a. Besides the volume 
of r is ^(a/n)'^, therefore 



I 



a ^, ,2 



|VP„(/U„)(a;)| dx = ^-/j l'"".^ ~ 1^' 



n,x\ 



For a center tetrahedron r* e T^, it is a bit more tricky since the edges are 
not orthogonal. The volume is of course the complement of the volumes of 
the corner tetrahedra, namely ^{a/n)^. The computation yields 



^2 



I * 



and hence 



I 



|VP„(/X„)(x)|2dx = 2 \^n,xf - t^n,xf I 



|2 
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I 



Gathering all the contributions 



reTn 



= 2 [ |VP„(^„)(X)|MX+ J] f \VPn{fin){x)\Mx 

4 ^ — l/ijixi l-^n,Xi\ ~l~ 2 y — - |/^n,Xi f^n,Xi\ '^n,Q- 

..-n-^^ 6ra ^.^^ 12n 

The coefficient 4 in front of the sum on En occurs because each element of 
En is an element of four t e Tn. Similarly each element of Cn belongs to two 
r* ^ T*. The positive error Sn,n is due to an over-estimation because some 
of the edges e e En are on the outer surface of fi„ and have been counted 
too many times. Following Remark HI the contribution of tS„ ^ will always 
be 0{l/n) less than that of the other terms, and therefore will tend to zero 
as n ^ cx). 

We rewrite the first sum 



A I |2 V~l (l I 

{i,j)eEn xeCnsi l/eAf„,n,j; 



(here each edge is counted twice through the couples {x,y) and (2/,a;)). For 
X G Zn,Q, we can only estimate {a/3n)\fin,x — /^n,j/P ^ 4a/3n, but there are 
only 0{n) such terms. For x e Cn,n\ln,n, (a/3'^)|/^n,x — /^n,yP ^ ac/3n^ and 
there are O(n^) such terms. Therefore 

4 S ^l/^n,x, -/^„,x/ =0(1). 

For the second sum 

V —\ — |2 — 1 V -^1 _ 4_ _ |2 

{i,j)eC„ {i,j,k)eS„ 



y —( 

, ^ 12nV 

(ij,fc)eS„ 



\^n,Xi l-^n,Xk\ ' \f^n,Xj l^n,Xk 



n,Xi f^n,Xk 
{i,j)eE„ (i,j,k)eS„ 
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As for the previous sum, we can decompose these sums into two 0(1) con- 
tributions. There are at most 0{n) terms contributing to Sn,n and stem- 
ming from an X e ln,n, therefore Sn,n = 0(1). Therefore ||Pn(^n)||Hi(n;R3) 
is uniformly bounded and Pn(/^n) is weakly convergent in H^(Q;M.^). This, 
together with our convergence definition [21 leads to Lemma [2l 

A by-product of this proof is the fact that we can write 

l|VPn(^n)||i2(f^.£(K3.K3)) = |VP„(/i„)(x)|Mx = ^£n,ex(^n) + Onif^n) 

where 

{i,j,k)eS„ 



3.3.2 Proof of Lemma [3t approximation of a function 

Following Ziemer theorem ([lO], Theorem 3.11.6), we know that for any func- 
tion u e if^(f2;]R^) and for all e > 0, there exists a function e C^(fi;]R^) 
such that X{uJu,e) ^ where 

Uu,e '■= {x e Q such that u{x) ^ Ue{x)}. 

We want to extend this result and be able to localize the irregularities of 
u with respect to a shrinking lattice. 

Let u e H^(Q; R^), and be the set of points where u is not C^. Since 
the gradient of m g /7^(r2;M^) has to be defined almost everywhere, there 

o 

cannot be an open ball in and therefore X^ = 0- 

We now fix e > 0. The Lebesgue and "H^ Hausdorff measures coincide in 
and therefore we both have X{Xu) ^ £ and ^^{Xu) ^ £■ 

Hence there exists a sequence of open balls {Bi)ie^ such that 

00 

Xu^l^Bi and ^ diam(Sj)^ < e. 

Since Xu is closed and bounded, it is compact and we can extract from 
this open cover a finite subcover. 

Hence there exists N e N such that 

TV TV 

Xu^l^Bi and ^diam(i3j)^ < e. 

i=l i=l 
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Obviously [Jj^g bounded with a piecewise boundary, hence 

n^co #(i:n,n) J 

Besides u^^n_^q_ e H^{Kjf^iBi;M?), and we therefore can choose Us such 
that uju,e ^f^iBi. Finally 

lim — ,, ' — < e. 



In our proof, we set Vie = fi\a;moo,e and begin to work on the restricted 
shrunk lattice Cn,ne = ^n,Qf]^e- We also denote Vn^n^ the subset of ele- 
ments X e Cn,ne such that #(A^n,n,x ^e) ^ 6, that is the set of nodes which 
are too close to dfls to have their 6 nearest neighbors in fi^. 

Since dfls is piecewise for all 5 > 0, we know that #Pn,n. = O^n^). 

3.3.3 Proof of Lemma |4t Construction 

Let X e C„ Q and y e Nnn,x- Since u e C"'^(f2; M^), Vm is bounded by some 
constant C on f2 and 



\Pn{u)x - Pn{u)y\'' ^ C^\x-y\'' ^ 



a\2 C^a^ 



Therefore Pn{u) e . Clearly if \u\ = 1 on f2, forall x e Cn,Q, \Pn{u)x\ = 1 
and Pn{u) e W^. 

Lemma [2] implies that there exists lioo £ H^iVt^M?) such that Pn{u) > 

Uoo- Since u e C^{fl;W^), Pn{pn{u)) converges towards u. This convergence 
is pointwise and even uniform on Q. Hence Uoo = u. 
Last 



^n,ex(Pn(M)) = - Yj Yj MPn{u)y - Pn{l 



a 

^ 2j 2j A\pn[u)y - Pn{U] 

^ Y Y A\u{y) - u{x)\' 



|2 



Since u e C^(Q; M^), for all x e Cn,n and all y e Nn^n,x 
\u{y) - u{x)\^ 



\y-x 



2 



V7 / ^ 

vu[x) 



\y - x\ 



2 

~l~ ^n,x,yi 
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and therefore for all x e Cn,n\Dn,n 

2 Hy) - u{x)\' = 2\Vu{x)\' + e,,,,, 

and the errors en,x are uniformly o(l) as n ^ go. Hence 

£n,eAPn{u)) = 2A V (-) ' ( | V^Z (x) ^ + e„,,) + O (-) , 

where the 0{l/n) stems from x e Vn^n- Since lini„^oo Xix£/:„ n^'^/'^)'^^"'^ ^ 
and the border of Vt is piecewise 

lim y f-V |Vm(x)|2 = r |VM(x)|2dx 

n^co ^ \nJ Jo 

and 

lim £n,e^{Pn{u)) = £oo,cy,{u). 
n^oo 

3.3.4 Proof of Lemma O Lower semi-continuity 

With Definition [2] for the convergence, Pnif^n) in H^{VL\M?), and 

lim^inf \Pn{l^n)\m{n;M?) ^ II"^oo||h1(Q;R3) 

(see [5], Proposition 2.3). 

Since we have assumed that VL is compact and has a piecewise bound- 
ary, we have already seen that the convergence is strong in L^(fi;]R^) and 
therefore 

Let us first fix e and work in fi^. Thanks to Proposition HJ we know that 
TTT-oo is continuous and 

This implies that 

li^i^^f (ll^-Pn(A''n)llL2(n^;£(M3.]Kd)) — || VP„ (p„ (mco) ) ||^2(j7^ .£(]R3.]Kd))^ ^ 0. 

According to Lemma [21 we can write 

||VP„(^n)||L2(Q_^.£(IR3.]jd)) = ^(^^,ex(^n) +"n(/^n))- 

18 



We therefore know that 

Im^f (^^^n,ex(^n) " ^^n,ex(Pn ("^oo)) + " <(p„(moo))^ ^ 0. 

We spht this into 5 parts, namely 

^^n,ex(/^n) " ^^n,ex(Pn("^oo)) = On + &n, 

"n(Mn) - a^(Pn("^oo)) = C„ + (i„ + e„ + /„, 

which can be expressed using tpn^x '■= fJ^n,x — moo{x), defined for all n e N* 
and X e Cn,n,- 

On = ^ ^ ^ li^n^x — 4'n,y\ i 

Cn = ^ 2] ("^oo(a^i) - "^co(a;fe)) • (^n,Xfe - ^n,a;,), 

^" ^ ^ 2 ('^ri,^, - '0n,x-J • (mooCxfc) - moo(a:j)), 

(ij,fc)eS„ 
(i,j,k)eSn 

fn = <Sn,nAPn{'moo)) - Sn,nAP'n)- 

We show below that e„ ^ a^, and 6^, c„, dn and tend to zero. This implies 
that 

liminf(a„(^„) - an{Pn{moo))) ^ -^liminf(f^ (/i„) - {pn{moo))) 

n^ao ZA n— >oo 

and therefore 

liminf(£^ (^„) - {pn{m^))) ^ 0. 

n—>(X> 

Lemma H implies that lim„^oo ^^,ex(Pn("^oo)) = ^ot,cx("^oo), hence 

liminf £^ (^„) ^ <?^,cx("^oo)- 

n^oo 

which ends the proof. 
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Proof of e„ ^ a„. For all Xi, Xj, Xk e vCn.n^, 

Since each couple (^, j) e is an element of 4 triples in S'„, we have 

(i,j,k)eSn (i,j)eE„ 

which is a much stronger result than e„ ^ a„. 

Proof of c„ and dn 0. Let {il'n,x — i^n,y) ■ (^oo(l/) — ^a^iz)) be one term 
of the sum in Cn and sei v = z — y. Then 

and in the same sum there is also a term (■0n,x ~ i^n,y) ■ {j^aa{y) — i^aaiy — 
except for y e Vn^n,. 

Now since ^ W^, there exists a constant ^ such that 

|'0n,j; ~ '071,2/ 1 ^ ^■!/'~) 
77/ 

except for y e /„,f7-, but since there are 0{n) such nodes, their contribution 
in Cn tends to 0. We also have 

Tl 

-{{rriaaiy) -'maaiy- ^)) + i^oaiv) " "T-oo (?/ + v))) 

a 

V 

> (Vmoo(y) - Vmoo(2/)) • -p-r = 0, 

n—>oa \v\ 

and therefore 

CL 111 
7- 2 {i^n,Xi-i^n,x^)-{mao{Xk)-mao{Xj)) = 0{-)0{n^)0{-)o{-) = o(l). 

Remark 5. When y e T>n,Q,^, we can only say that \'il)n,x — ^n,y\ ^ C^a/n 
(except on ln,nj ^^'^ |^cx)(?/) — "^oo(-2)| ^ Cmd/n, and since H^Vn^o,^ = 0{n'^), 
the contribution of these nodes in Cn is a 

0{n')Oi-)0{-)0{-) + 0(n)0(-)0(-) = O(-). 
nun n n n 

The sum dn is treated in the same way. 
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Proof of /„ 0. The quantity is a sum of 0{n^) terms reading 

like {a/6n)\tj.n,x, - /^n,x, P for (ij) e E^, or (a/12n)|/i„,^^ - iJ,n,xj\'^ for e 
Cn- By Hypothesis ini only 0(n) among these terms can be only o(l/n) and 
the others are 0{l/n^). Hence «Sn,r2g(//n) = o(l). 

On the other hand, the fact that Uoo e C^(f2e,M^) ensures that all the terms 
in Sn,n,{Pn{moo)) are 0(l/n^), and therefore is Sn,n^{pn{moo)) = 0{l/n). 

Proof of bn 0. We use the fact that C^in,; R^) is dense in H\n,; 
for the II • ||Hi(ne:;R3) norm. Let 77 > 0, there exists g such 
that II moo — '^7?||Hi(Qe;M3) ^ V- Then bn can be split into two contributions 
bn = b1 + b'^ where 

^" = ^ 1j 1j ~ ^v){x) - {moo - m^){y)) ■ {^n,x - i^n,y)- 

For the first term, we notice that 

= ^ I] 2 (rnr^ix) - mr,{y)) ■ 'llJn,x 
yeC„,si^ xeNn,ft,y 

As in the previous proof, we write y as x + v and 

K = - X; iimj^(x) - m^ix + v)) - {nij^ix - v) - mrj{x))) ■ ilJn,x- 

We estimate 

{mrj{x) — mri{x + v)) — (m.^(x — v) — mr^{x)) 

.0 pi 

= V mr,{x + tv) ■ V &t — V mr^{x + tv) ■ V dt 
Jt=l Jt=o 

2 

= -VmJx) ■ — + Vm„(x) • — + 0(\v\'^) ^ — , 
\v\ \v\ 

where the constant c only depends on the second derivative of m^, which is 
bounded on VL^. Thus 
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Thanks to the compact injection of iy^(Jl;M^) in L^(J1;R^), we deduce 
that this term vanishes. We obviously have 

\^n\^-^ X! 2 lirn^o - m^){x) - {m^ - m^){y)\. 
For any function u e C^{fl; M^), and y = x + v, if [x, y] a fl, then 



u{y) - u{x) 



Wu{x + tv) ■ V dt 

Vu{x) v + I (Vu{x + tv) - Vu{x)) ■ V dt 
Jt=i 



dt 



Applying this to nioa and m^, 
(moo - rnr,){y) - (ma, - m^)(x) = (Vmoo(x) - Vm^(x)) ■ v 

+ (Vmaoix + tv) — Vniaoix)) ■ V 

— {Vmr,{x + tv) — 'Vmr,{x)) ■ V dt. 
Jt=l 

★ By definition of m^, \\moa—'iTi^\\m(n-R3) ^ ?7, which implies that \\Vmoa{x) 
Vm^(x)||i2(f^^.]K3) ^ ry. Hence 



2 J] ( Vmoo (x) - Vm^ {x)).{y-x) 



1/2 



1/2 



Since 

^— 1 / \ I 

hm > — Vmoo(x)— Vm„( 



X 



\Vmaoix) 



Vm^(a;) l^dx, 



lim ^^-^ V V (Vmao{x) -Vmrj{x)).{y - x) 



n^cc n 



^ C\\Vmoo{x) - Vm^(a;) ||L2(n,;R3) ^ Crj. 
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* Let us treat the second contribution involving the function m^. Since 

rO rO rt 

{Vmr^{x + tv) -Vmr^{x)) ■ V dt = W^mn{x + t'v) ■ v <g)v dt'dt, 

Jt=l Jt=l Jt'=0 

and V^m„ is uniformly bounded on fig, we can estimate 



r° 1 

2 2 {\/mr,{x + t{y-x))-Vm^{x))-{y-x) dt = 0{-). 



★ In the last contribution, function moo is only C^. 



{Vniaoix + t{y - x)) -Vniaoix)) ■ {y - x) dt 

r „ ,,c^^ „ •Jt=\ 

i£ ,x 

2 Xi |(Vmoo(x + tt;) - Vmoo(x)) • i;| dt 

^ Y |(Vmoo(x + tt;)-Vmoo(x))-— I dt. 



< ^ 



bmce moo G M^), the elements of the sum converge uniformly towards 

as r?, ^ GO, and therefore the integral is o(n^). 

Hence lim„^oo l^'^l ^ and therefore lim„^oo |^ra| ^ ^5 for all > 0. This 
leads to lim„^oo l^nl =0, and the lemma is proved. 

3.3.5 Proof of Lemma [6] 

The terms that occur in £^„,ex(/^n) but not in £n,exif^n) are those involing 
couples {x,y) where one at least of the nodes belong to There are 

eO{n^) such nodes. Hence, by Hypothesis [3], the contribution of ln,n in 
^n,ex(/Jn) tcuds to zcro, and there exists c > such that for all {/i„}„gN* £ W^, 

lim |f^ex(/in) - ^n,ex(^n)| ^ C£ . 
n— >oo ' 

We have already estimated Pn{l^) on Vt\ftn,ei where Vtri,e is a polyhedral 
subset of VL which, like oj^ has an 0{e) Lebesgue measure. We already know 
that Pn(/Wn)|Q\Q„,, e H'^(n\nn,^\^^) and P„(m„)|^^ — m,oop\Q„, in the sense 
of {n\nn,e ; K^) , therefore 

|Vmoo(a;)|^ dx ^ liminf |VP„(/i„)(2;)|Mx ^ Ce. 
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liminf £:„ ex(/^n) = liminf(£:^ ^^(^„) + (^„,cx(/in) - ^LeAf^n))) 

n—>oa n— >oo 

^ lim inf cx(^n)) + 1™ {£n,cAP'n) - ex(A*n)) 
n— >0C' ' n— >0C' ' 

(moo) + liminf (£„,ex(/Un) - <?n,ex(/"n)) 

+ liminf(£:„,ex(/^n) - ^n,exM) 
^ ^oo,ex(moo) - (C + Ci)e. 

Since this holds for all £ > 0, we finally deduce that 

liminf 6:„,ex(A<n) ^ ^oo,ex(moo). 
n— >cx) 

4 Other energy contributions 

4.1 Magnet ost at ics: demagnetizing energy 

We can define a mapping hd : -L^(R''^; M"^) i-^ L^(]R"^; M'"^) by: for all u e 
L^(M^;M^), hd{u) is solution in the sense of distributions to 

rot /id (w) = 0, 
div hd(u) = — div?i. 

When u is the magnetization, hd{u) is the demagnetizing field. Its energy is 

^diu) = Y||/id(li)||L2(R3.M3). 

For u e L^(r2;M^), we denote by u the L^(M^;]R^) function which equals u 
inside and outside Q. Hence for a spin distribution /in e (S^)"^"-", the 
demagnetizing energy is defined by 



4.2 Total energy 

We define the total energy summing up the exchange and demagnetizing 
energies both in the lattice context: 

and the limit continuous one: 
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Theorem 2. Let {Hn)nem ^ . In the sense of the topology defined by 
Definition \E 



n— >oo 



Moreover, this convergence is compatible with the unit norm constraint. 
Proof. We have aheady shown that 

n— >oo 

Moreover : L^(M^;M^) L^(R^;M^) is hnear and continuous. We 
therefore choose {fin)nem ^ such that > e H^{Q]M.'^). This 

n— >oo 

means that the sequence (P„(/i„)) is weakly convergent in M'^) to /i. 

Hence 

PnifJ-n) > fJ- and P„{fi„) > II. 

n— >oo n— >oo 

In particular £'„,d(/Jn) ^ ^^oo,d(/i)- 

n— >oo 

Besides if g H^{Vl] M^), we know that J9n(/i) ^ in H^{Vt] M^), which 

n— >oo 

ends the proof of 

and Theorem [21 □ 



5 Conclusion 

In this paper, we prove a F-convergence result from a discrete description 
of ferromagnetic materials at the microscopic scale to the continuous one. 
This result has been shown thanks to a rigidity hypothesis on the lattice of 
magnetic moments. This modeling hypothesis is based on the Heisenberg 
interaction phenomenon and could be justified by a time multi-scale study. 
The new hypothesis would take into account the speed of the Heisenberg 
relaxation compared to the Larmor precession process. 

The results in this paper are the seed in order to address the micro- 
mesoscopic limit for dynamic processes to be able to better understand the 
dissipation phenomena involved in the mesoscopic Landau-Lifchitz system. 
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